Development and utility of a clinical research informatics application for participant recruitment and workflow management for a return of results pilot trial in familial hypercholesterolemia in the Million Veteran Program.

Objective: The development of clinical research informatics tools and workflow processes associated with re-engaging biobank participants has become necessary as genomic repositories increasingly consider the return of actionable research results. Materials and Methods: Here we describe the development and utility of an informatics application for participant recruitment and enrollment management for the Veterans Affairs Million Veteran Program Return Of Actionable Results Study, a randomized controlled pilot trial returning individual genetic results associated with familial hypercholesterolemia. Results: The application is developed in Python-Flask and was placed into production in November 2021. The application includes modules for chart review, medication reconciliation, participant contact and biospecimen logging, survey recording, randomization, and documentation of genetic counseling and result disclosure. Three primary users, a genetic counselor and two research coordinators, and 326 Veteran participants have been integrated into the system as of February 23, 2023. The application has successfully handled 3367 task requests involving greater than 95 000 structured data points. Specifically, application users have recorded 326 chart reviews, 867 recruitment telephone calls, 158 telephone-based surveys, and 61 return of results genetic counseling sessions, among other available study tasks. Conclusion: The development of usable, customizable, and secure informatics tools will become increasingly important as large genomic repositories begin to return research results at scale. Our work provides a proof-of-concept for developing and using such tools to aid in managing the return of results process within a national biobank.

Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization.

BACKGROUND: Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (CHD; PRSCHD) for 5 genetic ancestry groups. METHODS: We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding and continuous shrinkage priors (polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods) applied to summary statistics from the largest multi-ancestry genome-wide association study meta-analysis for CHD to date, including 1.1 million participants from 5 major genetic ancestry groups. Following training and optimization in the Million Veteran Program, we evaluated the best-performing PRSCHD in 176 988 individuals across 9 diverse cohorts. RESULTS: Multi-ancestry polygenic risk score for CHD developed using pruning and thresholding methods and polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods outperformed ancestry-specific Polygenic risk score for CHD developed using pruning and thresholding methods and polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods across a range of tuning values. Two best-performing multi-ancestry PRSCHD (ie, polygenic risk score for CHD developed using pruning and thresholding methods optimized using a multi-ancestry population and polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods optimized using a multi-ancestry population) and 1 ancestry-specific (PRSCSxEUR) were taken forward for validation. Polygenic risk score for CHD developed using pruning and thresholding methods (PT) optimized using a multi-ancestry population demonstrated the strongest association with CHD in individuals of South Asian genetic ancestry and European genetic ancestry (odds ratio per 1 SD [95% CI, 2.75 [2.41-3.14], 1.65 [1.59-1.72]), followed by East Asian genetic ancestry (1.56 [1.50-1.61]), Hispanic/Latino genetic ancestry (1.38 [1.24-1.54]), and African genetic ancestry (1.16 [1.11-1.21]). Polygenic risk score for CHD developed using ancestry-based continuous shrinkage methods optimized using a multi-ancestry population showed the strongest associations in South Asian genetic ancestry (2.67 [2.38-3.00]) and European genetic ancestry (1.65 [1.59-1.71]), lower in East Asian genetic ancestry (1.59 [1.54-1.64]), Hispanic/Latino genetic ancestry (1.51 [1.35-1.69]), and the lowest in African genetic ancestry (1.20 [1.15-1.26]). CONCLUSIONS: The use of summary statistics from a large multi-ancestry genome-wide meta-analysis improved the performance of PRSCHD in most ancestry groups compared with single-ancestry methods. Despite the use of one of the largest and most diverse sets of training and validation cohorts to date, improvement of predictive performance was limited in African genetic ancestry. This highlights the need for larger Genome-wide association study datasets of underrepresented populations to enhance the performance of PRSCHD.

A functional genomic framework to elucidate novel causal non-alcoholic fatty liver disease genes.

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver pathology in western countries, with serious public health consequences. Efforts to identify causal genes for NAFLD have been hampered by the relative paucity of human data from gold-standard magnetic resonance quantification of hepatic fat. To overcome insufficient sample size, genome-wide association studies using NAFLD surrogate phenotypes have been used, but only a small number of loci have been identified to date. In this study, we combined GWAS of NAFLD composite surrogate phenotypes with genetic colocalization studies followed by functional in vitro screens to identify bona fide causal genes for NAFLD. Approach & Results: We used the UK Biobank to explore the associations of our novel NAFLD score, and genetic colocalization to prioritize putative causal genes for in vitro validation. We created a functional genomic framework to study NAFLD genes in vitro using CRISPRi. Our data identify VKORC1, TNKS, LYPLAL1 and GPAM as regulators of lipid accumulation in hepatocytes and suggest the involvement of VKORC1 in the lipid storage related to the development of NAFLD. Conclusions: Complementary genetic and genomic approaches are useful for the identification of NAFLD genes. Our data supports VKORC1 as a bona fide NAFLD gene. We have established a functional genomic framework to study at scale putative novel NAFLD genes from human genetic association studies.

Diet Quality and Epigenetic Aging in the Women's Health Initiative.

BACKGROUND: Higher diet quality scores are associated with a lower risk for many chronic diseases and all-cause mortality; however, it is unclear if diet quality is associated with aging biology. OBJECTIVE: This study aimed to examine the association between diet quality and a measure of biological aging known as epigenetic aging. DESIGN: A cross-sectional data analysis was used to examine the association between three diet quality scores based on self-reported food frequency questionnaire data and five measures of epigenetic aging based on DNA methylation (DNAm) data from peripheral blood. PARTICIPANTS/SETTING: This study included 4,500 postmenopausal women recruited from multiple sites across the United States (1993-98), aged 50 to 79 years, with food frequency questionnaire and DNAm data available from the Women's Health Initiative baseline visit. MAIN OUTCOME MEASURES: Five established epigenetic aging measures were generated from HumanMethylation450 Beadchip DNAm data, including AgeAccelHannum, AgeAccelHorvath, AgeAccelPheno, AgeAccelGrim, and DunedinPACE. STATISTICAL ANALYSES PERFORMED: Linear mixed models were used to test for associations between three diet quality scores (Healthy Eating Index, Dietary Approaches to Stop Hypertension, and alternate Mediterranean diet scores) and epigenetic aging measures, adjusted for age, race and ethnicity, education, tobacco smoking, physical activity, Women's Health Initiative substudy from which DNAm data were obtained, and DNAm-based estimates of leukocyte proportions. RESULTS: Healthy Eating Index, Dietary Approaches to Stop Hypertension, and alternate Mediterranean diet scores were all inversely associated with AgeAccelPheno, AgeAccelGrim, and DunedinPACE (P < 0.05), with the largest effects with DunedinPACE. A one standard deviation increment in diet quality scores was associated with a decrement (ß ± SE) in DunedinPACE z score of -0.097 ± 0.014 (P = 9.70 x 10-13) for Healthy Eating Index, -0.107 ± 0.014 (P = 1.53 x 10-14) for Dietary Approaches to Stop Hypertension, and -0.068 ± 0.013 (P = 2.31 x 10-07) for the alternate Mediterranean diet. CONCLUSIONS: In postmenopausal women, diet quality scores were inversely associated with DNAm-based measures of biological aging, particularly DunedinPACE.

Mendelian randomization analyses suggest a causal role for circulating GIP and IL-1RA levels in homeostatic model assessment-derived measures of ß-cell function and insulin sensitivity in Africans without type 2 diabetes.

BACKGROUND: In vitro and in vivo studies have shown that certain cytokines and hormones may play a role in the development and progression of type 2 diabetes (T2D). However, studies on their role in T2D in humans are scarce. We evaluated associations between 11 circulating cytokines and hormones with T2D among a population of sub-Saharan Africans and tested for causal relationships using Mendelian randomization (MR) analyses. METHODS: We used logistic regression analysis adjusted for age, sex, body mass index, and recruitment country to regress levels of 11 cytokines and hormones (adipsin, leptin, visfatin, PAI-1, GIP, GLP-1, ghrelin, resistin, IL-6, IL-10, IL-1RA) on T2D among Ghanaians, Nigerians, and Kenyans from the Africa America Diabetes Mellitus study including 2276 individuals with T2D and 2790 non-T2D individuals. Similar linear regression models were fitted with homeostatic modelling assessments of insulin sensitivity (HOMA-S) and ß-cell function (HOMA-B) as dependent variables among non-T2D individuals (n = 2790). We used 35 genetic variants previously associated with at least one of these 11 cytokines and hormones among non-T2D individuals as instrumental variables in univariable and multivariable MR analyses. Statistical significance was set at 0.0045 (0.05/11 cytokines and hormones). RESULTS: Circulating GIP and IL-1RA levels were associated with T2D. Nine of the 11 cytokines and hormones (exceptions GLP-1 and IL-6) were associated with HOMA-S, HOMA-B, or both among non-T2D individuals. Two-stage least squares MR analysis provided evidence for a causal effect of GIP and IL-RA on HOMA-S and HOMA-B in multivariable analyses (GIP ~ HOMA-S ß = - 0.67, P-value = 1.88 × 10-6 and HOMA-B ß = 0.59, P-value = 1.88 × 10-5; IL-1RA ~ HOMA-S ß = - 0.51, P-value = 8.49 × 10-5 and HOMA-B ß = 0.48, P-value = 5.71 × 10-4). IL-RA was partly mediated via BMI (30-34%), but GIP was not. Inverse variance weighted MR analysis provided evidence for a causal effect of adipsin on T2D (multivariable OR = 1.83, P-value = 9.79 × 10-6), though these associations were not consistent in all sensitivity analyses. CONCLUSIONS: The findings of this comprehensive MR analysis indicate that circulating GIP and IL-1RA levels are causal for reduced insulin sensitivity and increased ß-cell function. GIP's effect being independent of BMI suggests that circulating levels of GIP could be a promising early biomarker for T2D risk. Our MR analyses do not provide conclusive evidence for a causal role of other circulating cytokines in T2D among sub-Saharan Africans.

Cardiorespiratory fitness and risk of heart failure with preserved ejection fraction.

AIMS: Preventive strategies for heart failure with preserved ejection fraction (HFpEF) include pharmacotherapies and lifestyle modifications. However, the association between cardiorespiratory fitness (CRF) assessed objectively by a standardized exercise treadmill test (ETT) and the risk of HFpEF has not been evaluated. Thus, we evaluated the association between CRF and HFpEF incidence. METHODS AND RESULTS: We assessed CRF in US Veterans (624 551 men; mean age 61.2 ± 9.7 years and 43 179 women; mean age 55.0 ± 8.9 years) by a standardized ETT performed between 1999 and 2020 across US Veterans Affairs Medical Centers. All had no evidence of heart failure or myocardial infarction prior to completion of the ETT. We assigned participants to one of five age- and gender-specific CRF categories (quintiles) based on peak metabolic equivalents (METs) achieved during the ETT and four categories based on CRF changes in those with two ETT evaluations (n = 139 434) ≥1.0 year apart. During a median follow-up of 10.1 years (interquartile range 6.0-14.3 years), providing 6 879 229 person-years, there were 16 493 HFpEF events with an average annual rate of 2.4 events per 1000 person-years. The adjusted risk of HFpEF decreased across CRF categories as CRF increased, independent of comorbidities. For fit individuals (≥10.5 METs) the hazard ratio (HR) was 0.48 (95% confidence interval [CI] 0.46-0.51) compared with least fit (≤4.9 METs; referent). Being unfit carried the highest risk (HR 2.88, 95% CI 2.67-3.11) of any other comorbidity. The risk of unfit individuals who became fit was 37% lower (HR 0.63, 95% CI 0.57-0.71), compared to those who remained unfit. CONCLUSIONS: Higher CRF levels are independently associated with lower HFpEF in a dose-response manner. Changes in CRF reflected proportional changes in HFpEF risk, suggesting that the HFpEF risk was modulated by CRF.

CYP2C19 Polymorphisms and Clinical Outcomes Following Percutaneous Coronary Intervention (PCI) in the Million Veterans Program.

Background: CYP2C19 loss-of-function (LOF) alleles decrease the antiplatelet effect of clopidogrel following percutaneous coronary intervention (PCI) in patients presenting with acute coronary syndrome (ACS). The impact of genotype in stable ischemic heart disease (SIHD) is unclear. Objectives: Determine the association of CYP2C19 genotype with major adverse cardiac events (MACE) after PCI for ACS or SIHD. Methods: Million Veterans Program (MVP) participants age <65 years with a PCI documented in the VA Clinical Assessment, Reporting and Tracking (CART) Program between 1/1/2009 to 9/30/2017, treated with clopidogrel were included. Time to MACE defined as the composite of all-cause death, stroke or myocardial infarction within 12 months following PCI. Results: Among 4,461 Veterans (mean age 59.1 ± 5.1 years, 18% Black); 44% had ACS, 56% had SIHD and 29% carried a CYP2C19 LOF allele. 301 patients (6.7%) experienced MACE while being treated with clopidogrel, 155 (7.9%) in the ACS group and 146 (5.9%) in the SIHD group. Overall, MACE was not significantly different between LOF carriers vs. noncarriers (adjusted hazard ratio [HR] 1.18, confidence interval [95%CI] 0.97-1.45, p=0.096). Among patients presenting with ACS, MACE risk in LOF carriers versus non-carriers was numerically higher (HR 1.30, 95%CI 0.98-1.73, p=0.067). There was no difference in MACE risk in patients with SIHD (HR 1.09, 95%CI 0.82-1.44; p=0.565). Conclusions: CYP2C19 LOF carriers presenting with ACS treated with clopidogrel following PCI experienced a numerically greater elevated risk of MACE events. CYP2C19 LOF genotype is not associated with MACE among patients presenting with SIHD.

CXCL12 regulates coronary artery dominance in diverse populations and links development to disease.

Mammalian cardiac muscle is supplied with blood by right and left coronary arteries that form branches covering both ventricles of the heart. Whether branches of the right or left coronary arteries wrap around to the inferior side of the left ventricle is variable in humans and termed right or left dominance. Coronary dominance is likely a heritable trait, but its genetic architecture has never been explored. Here, we present the first large-scale multi-ancestry genome-wide association study of dominance in 61,043 participants of the VA Million Veteran Program, including over 10,300 Africans and 4,400 Admixed Americans. Dominance was moderately heritable with ten loci reaching genome wide significance. The most significant mapped to the chemokine CXCL12 in both Europeans and Africans. Whole-organ imaging of human fetal hearts revealed that dominance is established during development in locations where CXCL12 is expressed. In mice, dominance involved the septal coronary artery, and its patterning was altered with Cxcl12 deficiency. Finally, we linked human dominance patterns with coronary artery disease through colocalization, genome-wide genetic correlation and Mendelian Randomization analyses. Together, our data supports CXCL12 as a primary determinant of coronary artery dominance in humans of diverse backgrounds and suggests that developmental patterning of arteries may influence one's susceptibility to ischemic heart disease.

A Multi-Ancestry Polygenic Risk Score for Coronary Heart Disease Based on an Ancestrally Diverse Genome-Wide Association Study and Population-Specific Optimization.

Background: Predictive performance of polygenic risk scores (PRS) varies across populations. To facilitate equitable clinical use, we developed PRS for coronary heart disease (PRSCHD) for 5 genetic ancestry groups. Methods: We derived ancestry-specific and multi-ancestry PRSCHD based on pruning and thresholding (PRSP+T) and continuous shrinkage priors (PRSCSx) applied on summary statistics from the largest multi-ancestry genome-wide meta-analysis for CHD to date, including 1.1 million participants from 5 continental populations. Following training and optimization of PRSCHD in the Million Veteran Program, we evaluated predictive performance of the best performing PRSCHD in 176,988 individuals across 9 cohorts of diverse genetic ancestry. Results: Multi-ancestry PRSP+T outperformed ancestry specific PRSP+T across a range of tuning values. In training stage, for all ancestry groups, PRSCSx performed better than PRSP+T and multi-ancestry PRS outperformed ancestry-specific PRS. In independent validation cohorts, the selected multi-ancestry PRSP+T demonstrated the strongest association with CHD in individuals of South Asian (SAS) and European (EUR) ancestry (OR per 1SD[95% CI]; 2.75[2.41-3.14], 1.65[1.59-1.72]), followed by East Asian (EAS) (1.56[1.50-1.61]), Hispanic/Latino (HIS) (1.38[1.24-1.54]), and weakest in African (AFR) ancestry (1.16[1.11-1.21]). The selected multi-ancestry PRSCSx showed stronger associacion with CHD in comparison within each ancestry group where the association was strongest in SAS (2.67[2.38-3.00]) and EUR (1.65[1.59-1.71]), progressively decreasing in EAS (1.59[1.54-1.64]), HIS (1.51[1.35-1.69]), and lowest in AFR (1.20[1.15-1.26]). Conclusions: Utilizing diverse summary statistics from a large multi-ancestry genome-wide meta-analysis led to improved performance of PRSCHD in most ancestry groups compared to single-ancestry methods. Improvement of predictive performance was limited, specifically in AFR and HIS, despite use of one of the largest and most diverse set of training and validation cohorts to date. This highlights the need for larger GWAS datasets of AFR and HIS individuals to enhance performance of PRSCHD.

A multi-ancestry polygenic risk score improves risk prediction for coronary artery disease.

Identification of individuals at highest risk of coronary artery disease (CAD)-ideally before onset-remains an important public health need. Prior studies have developed genome-wide polygenic scores to enable risk stratification, reflecting the substantial inherited component to CAD risk. Here we develop a new and significantly improved polygenic score for CAD, termed GPSMult, that incorporates genome-wide association data across five ancestries for CAD (>269,000 cases and >1,178,000 controls) and ten CAD risk factors. GPSMult strongly associated with prevalent CAD (odds ratio per standard deviation 2.14, 95% confidence interval 2.10-2.19, P < 0.001) in UK Biobank participants of European ancestry, identifying 20.0% of the population with 3-fold increased risk and conversely 13.9% with 3-fold decreased risk as compared with those in the middle quintile. GPSMult was also associated with incident CAD events (hazard ratio per standard deviation 1.73, 95% confidence interval 1.70-1.76, P < 0.001), identifying 3% of healthy individuals with risk of future CAD events equivalent to those with existing disease and significantly improving risk discrimination and reclassification. Across multiethnic, external validation datasets inclusive of 33,096, 124,467, 16,433 and 16,874 participants of African, European, Hispanic and South Asian ancestry, respectively, GPSMult demonstrated increased strength of associations across all ancestries and outperformed all available previously published CAD polygenic scores. These data contribute a new GPSMult for CAD to the field and provide a generalizable framework for how large-scale integration of genetic association data for CAD and related traits from diverse populations can meaningfully improve polygenic risk prediction.

Diversity and Scale: Genetic Architecture of 2,068 Traits in the VA Million Veteran Program.

Genome-wide association studies (GWAS) have underrepresented individuals from non-European populations, impeding progress in characterizing the genetic architecture and consequences of health and disease traits. To address this, we present a population-stratified phenome-wide GWAS followed by a multi-population meta-analysis for 2,068 traits derived from electronic health records of 635,969 participants in the Million Veteran Program (MVP), a longitudinal cohort study of diverse U.S. Veterans genetically similar to the respective African (121,177), Admixed American (59,048), East Asian (6,702), and European (449,042) superpopulations defined by the 1000 Genomes Project. We identified 38,270 independent variants associating with one or more traits at experiment-wide P<4.6×10-11 significance; fine-mapping 6,318 signals identified from 613 traits to single-variant resolution. Among these, a third (2,069) of the associations were found only among participants genetically similar to non-European reference populations, demonstrating the importance of expanding diversity in genetic studies. Our work provides a comprehensive atlas of phenome-wide genetic associations for future studies dissecting the architecture of complex traits in diverse populations.

Systems Age: A single blood methylation test to quantify aging heterogeneity across 11 physiological systems.

Individuals, organs, tissues, and cells age in diverse ways throughout the lifespan. Epigenetic clocks attempt to quantify differential aging between individuals, but they typically summarize aging as a single measure, ignoring within-person heterogeneity. Our aim was to develop novel systems-based methylation clocks that, when assessed in blood, capture aging in distinct physiological systems. We combined supervised and unsupervised machine learning methods to link DNA methylation, system-specific clinical chemistry and functional measures, and mortality risk. This yielded a panel of 11 system-specific scores- Heart, Lung, Kidney, Liver, Brain, Immune, Inflammatory, Blood, Musculoskeletal, Hormone, and Metabolic. Each system score predicted a wide variety of outcomes, aging phenotypes, and conditions specific to the respective system, and often did so more strongly than existing epigenetic clocks that report single global measures. We also combined the system scores into a composite Systems Age clock that is predictive of aging across physiological systems in an unbiased manner. Finally, we showed that the system scores clustered individuals into unique aging subtypes that had different patterns of age-related disease and decline. Overall, our biological systems based epigenetic framework captures aging in multiple physiological systems using a single blood draw and assay and may inform the development of more personalized clinical approaches for improving age-related quality of life.

Plasma proteomic signatures of a direct measure of insulin sensitivity in two population cohorts.

AIMS/HYPOTHESIS: The euglycaemic-hyperinsulinaemic clamp (EIC) is the reference standard for the measurement of whole-body insulin sensitivity but is laborious and expensive to perform. We aimed to assess the incremental value of high-throughput plasma proteomic profiling in developing signatures correlating with the M value derived from the EIC. METHODS: We measured 828 proteins in the fasting plasma of 966 participants from the Relationship between Insulin Sensitivity and Cardiovascular disease (RISC) study and 745 participants from the Uppsala Longitudinal Study of Adult Men (ULSAM) using a high-throughput proximity extension assay. We used the least absolute shrinkage and selection operator (LASSO) approach using clinical variables and protein measures as features. Models were tested within and across cohorts. Our primary model performance metric was the proportion of the M value variance explained (R2). RESULTS: A standard LASSO model incorporating 53 proteins in addition to routinely available clinical variables increased the M value R2 from 0.237 (95% CI 0.178, 0.303) to 0.456 (0.372, 0.536) in RISC. A similar pattern was observed in ULSAM, in which the M value R2 increased from 0.443 (0.360, 0.530) to 0.632 (0.569, 0.698) with the addition of 61 proteins. Models trained in one cohort and tested in the other also demonstrated significant improvements in R2 despite differences in baseline cohort characteristics and clamp methodology (RISC to ULSAM: 0.491 [0.433, 0.539] for 51 proteins; ULSAM to RISC: 0.369 [0.331, 0.416] for 67 proteins). A randomised LASSO and stability selection algorithm selected only two proteins per cohort (three unique proteins), which improved R2 but to a lesser degree than in standard LASSO models: 0.352 (0.266, 0.439) in RISC and 0.495 (0.404, 0.585) in ULSAM. Reductions in improvements of R2 with randomised LASSO and stability selection were less marked in cross-cohort analyses (RISC to ULSAM R2 0.444 [0.391, 0.497]; ULSAM to RISC R2 0.348 [0.300, 0.396]). Models of proteins alone were as effective as models that included both clinical variables and proteins using either standard or randomised LASSO. The single most consistently selected protein across all analyses and models was IGF-binding protein 2. CONCLUSIONS/INTERPRETATION: A plasma proteomic signature identified using a standard LASSO approach improves the cross-sectional estimation of the M value over routine clinical variables. However, a small subset of these proteins identified using a stability selection algorithm affords much of this improvement, especially when considering cross-cohort analyses. Our approach provides opportunities to improve the identification of insulin-resistant individuals at risk of insulin resistance-related adverse health consequences.

Cardiovascular Disease Risk Assessment Using Traditional Risk Factors and Polygenic Risk Scores in the Million Veteran Program.

Importance: Primary prevention of atherosclerotic cardiovascular disease (ASCVD) relies on risk stratification. Genome-wide polygenic risk scores (PRSs) are proposed to improve ASCVD risk estimation. Objective: To determine whether genome-wide PRSs for coronary artery disease (CAD) and acute ischemic stroke improve ASCVD risk estimation with traditional clinical risk factors in an ancestrally diverse midlife population. Design, Setting, and Participants: This was a prognostic analysis of incident events in a retrospectively defined longitudinal cohort conducted from January 1, 2011, to December 31, 2018. Included in the study were adults free of ASCVD and statin naive at baseline from the Million Veteran Program (MVP), a mega biobank with genetic, survey, and electronic health record data from a large US health care system. Data were analyzed from March 15, 2021, to January 5, 2023. Exposures: PRSs for CAD and ischemic stroke derived from cohorts of largely European descent and risk factors, including age, sex, systolic blood pressure, total cholesterol, high-density lipoprotein (HDL) cholesterol, smoking, and diabetes status. Main Outcomes and Measures: Incident nonfatal myocardial infarction (MI), ischemic stroke, ASCVD death, and composite ASCVD events. Results: A total of 79 151 participants (mean [SD] age, 57.8 [13.7] years; 68 503 male [86.5%]) were included in the study. The cohort included participants from the following harmonized genetic ancestry and race and ethnicity categories: 18 505 non-Hispanic Black (23.4%), 6785 Hispanic (8.6%), and 53 861 non-Hispanic White (68.0%) with a median (5th-95th percentile) follow-up of 4.3 (0.7-6.9) years. From 2011 to 2018, 3186 MIs (4.0%), 1933 ischemic strokes (2.4%), 867 ASCVD deaths (1.1%), and 5485 composite ASCVD events (6.9%) were observed. CAD PRS was associated with incident MI in non-Hispanic Black (hazard ratio [HR], 1.10; 95% CI, 1.02-1.19), Hispanic (HR, 1.26; 95% CI, 1.09-1.46), and non-Hispanic White (HR, 1.23; 95% CI, 1.18-1.29) participants. Stroke PRS was associated with incident stroke in non-Hispanic White participants (HR, 1.15; 95% CI, 1.08-1.21). A combined CAD plus stroke PRS was associated with ASCVD deaths among non-Hispanic Black (HR, 1.19; 95% CI, 1.03-1.17) and non-Hispanic (HR, 1.11; 95% CI, 1.03-1.21) participants. The combined PRS was also associated with composite ASCVD across all ancestry groups but greater among non-Hispanic White (HR, 1.20; 95% CI, 1.16-1.24) than non-Hispanic Black (HR, 1.11; 95% CI, 1.05-1.17) and Hispanic (HR, 1.12; 95% CI, 1.00-1.25) participants. Net reclassification improvement from adding PRS to a traditional risk model was modest for the intermediate risk group for composite CVD among men (5-year risk >3.75%, 0.38%; 95% CI, 0.07%-0.68%), among women, (6.79%; 95% CI, 3.01%-10.58%), for age older than 55 years (0.25%; 95% CI, 0.03%-0.47%), and for ages 40 to 55 years (1.61%; 95% CI, -0.07% to 3.30%). Conclusions and Relevance: Study results suggest that PRSs derived predominantly in European samples were statistically significantly associated with ASCVD in the multiancestry midlife and older-age MVP cohort. Overall, modest improvement in discrimination metrics were observed with addition of PRSs to traditional risk factors with greater magnitude in women and younger age groups.

Mendelian randomization analyses clarify the effects of height on cardiovascular diseases.

An inverse correlation between stature and risk of coronary artery disease (CAD) has been observed in several epidemiologic studies, and recent Mendelian randomization (MR) experiments have suggested causal association. However, the extent to which the effect estimated by MR can be explained by established cardiovascular risk factors is unclear, with a recent report suggesting that lung function traits could fully explain the height-CAD effect. To clarify this relationship, we utilized a well-powered set of genetic instruments for human stature, comprising >1,800 genetic variants for height and CAD. In univariable analysis, we confirmed that a one standard deviation decrease in height (~6.5 cm) was associated with a 12.0% increase in the risk of CAD, consistent with previous reports. In multivariable analysis accounting for effects from up to 12 established risk factors, we observed a >3-fold attenuation in the causal effect of height on CAD susceptibility (3.7%, p = 0.02). However, multivariable analyses demonstrated independent effects of height on other cardiovascular traits beyond CAD, consistent with epidemiologic associations and univariable MR experiments. In contrast with published reports, we observed minimal effects of lung function traits on CAD risk in our analyses, indicating that these traits are unlikely to explain the residual association between height and CAD risk. In sum, these results suggest the impact of height on CAD risk beyond previously established cardiovascular risk factors is minimal and not explained by lung function measures.

Epigenome-wide meta-analysis of BMI in nine cohorts: Examining the utility of epigenetically predicted BMI.

This study sought to examine the association between DNA methylation and body mass index (BMI) and the potential of BMI-associated cytosine-phosphate-guanine (CpG) sites to provide information about metabolic health. We pooled summary statistics from six trans-ethnic epigenome-wide association studies (EWASs) of BMI representing nine cohorts (n = 17,034), replicated these findings in the Women's Health Initiative (WHI, n = 4,822), and developed an epigenetic prediction score of BMI. In the pooled EWASs, 1,265 CpG sites were associated with BMI (p < 1E-7) and 1,238 replicated in the WHI (FDR < 0.05). We performed several stratified analyses to examine whether these associations differed between individuals of European and African descent, as defined by self-reported race/ethnicity. We found that five CpG sites had a significant interaction with BMI by race/ethnicity. To examine the utility of the significant CpG sites in predicting BMI, we used elastic net regression to predict log-normalized BMI in the WHI (80% training/20% testing). This model found that 397 sites could explain 32% of the variance in BMI in the WHI test set. Individuals whose methylome-predicted BMI overestimated their BMI (high epigenetic BMI) had significantly higher glucose and triglycerides and lower HDL cholesterol and LDL cholesterol compared to accurately predicted BMI. Individuals whose methylome-predicted BMI underestimated their BMI (low epigenetic BMI) had significantly higher HDL cholesterol and lower glucose and triglycerides. This study confirmed 553 and identified 685 CpG sites associated with BMI. Participants with high epigenetic BMI had poorer metabolic health, suggesting that the overestimation may be driven in part by cardiometabolic derangements characteristic of metabolic syndrome.

Genetic evidence for causal relationships between age at natural menopause and the risk of ageing-associated adverse health outcomes.

BACKGROUND: A later age at natural menopause (ANM) has been linked to several ageing-associated traits including an increased risk of breast and endometrial cancer and a decreased risk of lung cancer, osteoporosis and Alzheimer disease. However, ANM is also related to several proxies for overall health that may confound these associations. METHODS: We investigated the causal association of ANM with these clinical outcomes using Mendelian randomization (MR). Participants and outcomes analysed were restricted to post-menopausal females. We conducted a one-sample MR analysis in both the Women's Health Initiative and UK Biobank. We further analysed and integrated several additional data sets of post-menopausal women using a two-sample MR design. We used ≤55 genetic variants previously discovered to be associated with ANM as our instrumental variable. RESULTS: A 5-year increase in ANM was causally associated with a decreased risk of osteoporosis [odds ratio (OR) = 0.80, 95% CI (0.70-0.92)] and fractures (OR = 0.76, 95% CI, 0.62-0.94) as well as an increased risk of lung cancer (OR = 1.35, 95% CI, 1.06-1.71). Other associations including atherosclerosis-related outcomes were null. CONCLUSIONS: Our study confirms that the decline in bone density with menopause causally translates into fractures and osteoporosis. Additionally, this is the first causal epidemiological analysis to our knowledge to find an increased risk of lung cancer with increasing ANM. This finding is consistent with molecular and epidemiological studies suggesting oestrogen-dependent growth of lung tumours.

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program.

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3. Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.

Genomics and phenomics of body mass index reveals a complex disease network.

Elevated body mass index (BMI) is heritable and associated with many health conditions that impact morbidity and mortality. The study of the genetic association of BMI across a broad range of common disease conditions offers the opportunity to extend current knowledge regarding the breadth and depth of adiposity-related diseases. We identify 906 (364 novel) and 41 (6 novel) genome-wide significant loci for BMI among participants of European (N~1.1 million) and African (N~100,000) ancestry, respectively. Using a BMI genetic risk score including 2446 variants, 316 diagnoses are associated in the Million Veteran Program, with 96.5% showing increased risk. A co-morbidity network analysis reveals seven disease communities containing multiple interconnected diseases associated with BMI as well as extensive connections across communities. Mendelian randomization analysis confirms numerous phenotypes across a breadth of organ systems, including conditions of the circulatory (heart failure, ischemic heart disease, atrial fibrillation), genitourinary (chronic renal failure), respiratory (respiratory failure, asthma), musculoskeletal and dermatologic systems that are deeply interconnected within and across the disease communities. This work shows that the complex genetic architecture of BMI associates with a broad range of major health conditions, supporting the need for comprehensive approaches to prevent and treat obesity.